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Abstract 

We study the geodesic motion in Godel's universe, using conserved quantities. We give a 
necessary and suSicient condition for curves to be geodesic curves in terms of conserved quantities, 
which can be computed from the initial values of the curve. We check our result with numerical 
simulations too. 

1 Introduction 

Godel's solution [5] for the Einstein equation gives a cosmological model of a rotating universe. This 
solution has many interesting properties, for example it contains closed timelike curves, geodesically 
complete, and has neither a singularity nor a horizon. The free motion of particles in this cosmological 
model was analyzed first by Kundt [B] and Chandrasekhar and Wright [I]. A physical model for the 
geodesic motion was given byNovelloet. al. [7] using the effective potential. In this paper we describe 
the free motion with conserved quantities. 

In the flat spacctimc (M^, <? = Diag(l, —1, —1, —1)) the equation for a geodesic curve 7 : M ^ R'' 
is given by 

7 = 0. (1) 

The solution is obvious "/{u) — vt + x. The vectors u, a; are given as initial values. In this case we 
can consider vas a conserved quantity, since a curve 7 describes a free motion if and only if 'j{u) = v 
for every u G Dom(7). In this setting the physical meaning of v is well-known. In the next section 
we give a similar description for geodesic motion in the Godel spacetime and we check our result 
using numerical simulations. During the computation we use the Einstein summation and the usual 
differential geometrical formalism. 

2 Momentums in the Godel's universe 



2.1 Geodesic lines 

Godel presented the line-element of his universe in the following form [2], [31 P- 195], [H p. 275]. 
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Deflnition 2.1. We call the Riemannian manifold {Ai,g) to Gddel spacetime, where the manifold is 
A4 expressed in cylindrical coordinates M x M+ x [0, 27r[ x R and the Riemannian metric at a point 
{t, r, ip,z) e M is 



/ 1 %/2sinhV 0\ 

0-1 

\/2sinh^(r) sinh'*(r) - sinh^M 

V -1/ 



(2) 



where a is positive parameter. 

The parameter a can be interpreted as a = "^^^i where w is a measure of the constant rotation 
H p. 191]. 

To get the differential equations of the geodesic curve we compute the second order Christoffel 
symbols, using the Equation 

r™ = + d.g.k - dug,,) , (3) 

where g^™ denotes the {k,m) element of the inverse matrix of g. The nonzero symbols are the 

'k pfc 

ij ~ ji 



following (and their counterparts F*^ — F^^J 



, 2sinh(r) , _ v/2sinh^(r) 

cosh(r) ^23- ^Qgj^(^) W 



2 V2sinh(r) cosh(r) , sinh(r) cosh(r)(2 cosh2(r) — 3) 

^13= ^23= (5) 



sinh(r) cosh(r) sinh(r) cosh(r) 

A curve 7 : R ^ Af is a geodesic line if 



7'+r^,f7' = (7) 

holds for k — 1,2, 3, 4. If we write the curve as 7(u) = {t{u), r{u), <p{u), z{u)) then 7 is a geodesic line 
if and only if the following equations hold. 

.. 4sinh(r) . 2V2sinh^(r) , , 

t + -^tr + ; ' ^r = 8 

cosh(r) cosh(r) 

r + 2\/2sinh(r) cosh(r)i(^ + sinh(r) cosh(r)(2 cosh2(r) - 3)((^)2 = (9) 

V? sinh(r) cosh(r) - 2V2if + 2(^f = (10) 

2 = 0, (11) 

where the dot ' denotes the derivation with respect to u. It is known from the theory of differential 
equations that, there exists a unique solution for every initial value 0. 

2.2 Description of the free motion with conserved quantities 

In this subsection we derive a 4 x 4 matrix M, which mixes the components of 7, where 7 is a curve in 
the Godel spacetime, such that 7 is a geodesic curve, if and only if A/7 is constant. Since during the 
geodesic motion Af7 is constant, like the impulse momentum in the classical case, we call this matrix 
M to momentum matrix. First we define this matrix, and then we prove its mentioned properties. 
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Definition 2.2. For parameters pi,p2 G R, r G and (p E [0, 27r[ let us define the momentum 
matrix M(pi, p2, r, ip) 



I picosli^(r) 
Pi H 7, 



— \/2 cos(i^) sinh(2r) sin((^) 
\/2 sin((p) sinh(2r) eos((/7) 

V 



(pi cosh^ (r ) + 4p2 ) sinh^ (r ) \ 



2V2 




eos(iy9)(eosli(2r) - 


- 2) sinfi(2r) 


2 




sin((/9)(cosh(2r) — 


2) sinh(2r) 


2 
















(12) 



Theorem 2.3. for every geodesic line "f{u) — {t{u),r(u), if{u), z{u)) and for every real parameters 
Pi, P2 o-nd there exist real parameters Ci,C2,C3 and f{pi,P2) depending just on the initial conditions 
7(mo) and jiuo), such that 



M{pi,p2,r, (f) 



r 



ffiPUP2)\ 
Cl 
C2 

V C3 / 



(13) 



holds at every point of Donry. If a curve 7 : R ^ Af satisfies Equation US\) at every point o/Dom7, 
then 7 is a geodesic curve. 

Proof. Since the Theorem is obviously valid for the last component z, we skip it during the proof. 
Consider the left hand side of the differential Equations (0 El IIHl) as a vector X = (dHl), ®, (fTU| ). 
Assume that D{t^ r, (^) is a 3 x 3 is an invertible matrix, which entries are smooth functions of t, r and 
if. Assume that the equation 



d 

d u 



(^D{t{u),r{u),ipiu)) ■ 7(u)) = D{t{u), r{u), <f{u)) ■ X(ti) 



(14) 



holds. 



If D consists of the column- vectors di , ^2 and ^3 , then the Equation can be written as 

cht + d2r + dsp + -^{tf + -^{ff + -={^' + 

ddi dd2\ . fddi 9^3 \ • f dd2 dd3\ 



(15) 



dii' + d2r + fi3</3 + sinh(r) cosh(r)(2 cosh^(r) - 3)rf2(<^)^ + 

-ds if + 2\/2 sinh(r) cosh(r)(i2^u</5+ 



^ ^ sinh(r) 2V2 

4 — — rdi — — 



cosh(r ) — sinh(r) cosh(r) — ; 

2 



2V2'^^d, 



cosh(r) — sinh(r) cosh(r) — 



^3 Vr- 
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After elementary calculations wc have the set of differential equations for the vectors ((ii)i=i,2,3- 

(16) 



'dt 



9di dd2 2V2 , . , w ^ 

or at sinh(r) cosh(rj — — 

ddi dd3 2sinh(2r), , 

dip dt y/2 — 

(19) 
da + 2\/2cosh(r) sinh(r)di - 2\/2di tanh(r) (20) 



Or 

dd2 dds 



dtp dr sinh(r) cosh(r)) 
dd 1 

-= = -d2(cosh(2r) - 2) sinh(2r) (21) 

Let us note, that the same set of differential equations are valid for all components of the vectors 
di,d2,d3. Therefore it is enough to consider just the first component of the vectors, which will be 
denoted with di,d2 and d^. (That is for index i £ {1,2,3}, di denotes the first component of the 
vector di.) 

Equations p6ll9p imply that di = di(r, ip) and d2 ~ ^2(^1 'P')- If we compute the derivative of the 
Equation (jlSp with respect to ip and t, using Equation (j2ip we have the following equation for 1^2 . 

/ cosh(2r)-2 \ d^d2 r- dH2 

i 2 j^^^^dW^ (22^ 

Computing the derivative of the previous Equation with respect to r, and taking into account the 
Equation (fT9|) we have the following form for ^2- 

d2{t,ip)^bi{ip)t + b2{ip) (23) 

Using this form for d2 Equation (|2ip gives 

d3(t,r,v3) = i (i / bx{ip)Aip+ / h2{ip)Aip] •(cosh(2r)-2)sinh(2r)+e(r,t), (24) 



where e{r,t) is a smooth function. Substituting this form of ^3 to the Equation (|17p and computing 
the derivative with respect t, we have 

= -2 V2(cosh(2r ) - 2) / 5i (^) d ^ . (25) 

J cosh(r) smh(r) ot 

This means, that there exists a real number A € R, such that 

fbi{^)d^^X, and — — J ^^jfi^ ^ -A (26) 

y ^ (cosh(2r) - 2) cosh(r) sinh(r) dt ^ ' 
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holds. This means that bi{ifi) — 0. From this d2 = 62(9') and e{r,t) — g{r) foUows. So far, we have 
the foUowing forms for the functions ((ii)i=i,2,3- 

di= V2smh{2r) [ b2iip) dip + f{r) (27) 



d2 = b2{^) (28) 

= ^ (/ ^2((^) d(^^ (cosh(2r) - 2) sinh(2r) + g{r) (29) 

Where Equation ([27]) comes from Equation (fTS]) . Let us substitute the Equations P7l28l29p into the 
Equation 



^ = / 62(<^)d^^ f2(cosh(2r)-2)cosh(2r) + 2sinh2(2r)') -4^(r) 

aip 2 \J J \ / dr 



62 (t/') d (/3 (cosh(2r) - 2) sinh(2r) + g{r) 



cosh(r) sinh(r) \2 
+ 2V2 ( \/2sinh(2r) / 62(¥') d(y5 + /(r) ) cosh(r) sinh(r) 



-2^(2) f^/2sinh(2r) J b^i^) d ^ + f {r)^ 



After simphfications we get the foUowing equation 
d62((p) 



d If 

dg{r) , 2g{r 



d r cosh(r) sinh(r) 

This imphes that 



b2i(p)dip 

+ 2V2/(r) cosh(r) sinh(r) - 2%/2/(r) tanh(r) 



(30) 



(31) 



^^^^^)+ /62(^)d^ = a (32) 



dip 

holds, for a real constant a. This gives the explicit form of 62 

62(95) = ki cos((p) + k2 sin(93), (33) 
where fci, /c2 G K- After this form of 62 we rewrite the Equations (|27l28l29p . 

di = V2(^ki sm{ip) - k2 cos(^)) sinh(2r) + /*(r) (34) 

d2 — ki cos{ip) + k2 sin(iy9) (35) 

= i ^ki sm{<p) - k2 cos{ip)^ (^cosh(2r) - 2^ sinh(2r) + g*{r) (36) 
If we substitute Equations (|34|35l36p into Equation (|17p. after some simplification we have 

" • 1,(0 N + 4/ (r) tanh(r). (37) 
dr smh(zr) 



And if we substitute Equations (|34l35l36p into Equation (1^0]) . then after simplifications we get 

d<?% , 4(g*(r)+V2r(r)sinh^(r) 

If we express tfie function g*{r) from Equation (I57|) and substitute it to the Equation then we 
get a solvable differential equation. 

2 cosh(2r) ^ (r) = sinh(2r) (r) (39) 
dr dr"^ 

The general solution is 

r(r) =p2 + ^PiC0sh(rf, (40) 

where pi, P2 G After this, we can compute the function g* . 

(pi +8p2 + PiCosh(2r))sinh2(r) 
5 M = ^^^^ . (41) 

We determined the first components of the vectors di, d2 and d^. This means, that the column- vectors 
of the matrix D are of the following form. 

di ^ \/2(fci sin((^) - fc2 cos(v3)) sinh(2r) + l{p2 + ^Pi cosh^(r)) (42) 

^2 ~ ki cos{(p) + k2 sm((p) (43) 

ds = ^(fci sin((^) - A^cos(.^))(cosh(2r) - 2) sinh(2r) + (44) 

^ ^ (.91 + 8ff2 + 51 cosh(2r)) sinh^(r) 



4V2 

where 1 = ("TJ' "75' "73) ^^^^ ^^"^ vectors fci, fc2 and 1 are linearly independent, since we assumed, that 
the matrix D is invertible. Moreover we can assume that the vectors ki, k2 and 1_ form an orthonormal 
basis in R^. According to this, if we multiply from left the matrix D with the matrix L which has 
row-vectors (1, fe^, we get the matrix M{pi, p2,r,(p), which is independent of t. This proves the 
Theorem. □ 

2.3 Numerical verification of the Theorem 

Now we give a numerical verification of the previous Theorem. First we note that enough to consider 
two linearly independent constants vectors {pi,p2) & and (p'1,/32) ^ the matrix M, since 

easy to check that the function /(pi, P2) is linear 

f{Pl,P2)+.f{p'l,P2) = fiPl+Pl,P2+p'2)- (45) 

We choose the matrices M{0,l,t,r,(p) and M ^1, — -, t, r, (/3^ . For these M matrices the Equation 
()13p can be written in the following form. (We skip the z component, since this part of the Theorem 
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is obvious.) 



i(u) + V2sinh^(r)^(M) = /(0, 1) (46) 



cosh(2r) • sinh''(r) / 1 



4 ^(-)+^^(")-H''"iJ ^''^ 

— •\/2cos((p) sinh(2r)i(M) + sin{ip)f{u) 

cos((^)(cosh(2r)) ~ 2) sinh(2r) ^^^^^ ^ ^^^^ 

\/2ain{ip) sinh(2r)i(w) + cos(i^)r(ii) 

sin(y^)(cosh(2r) - 2)sinh(2r) 
+ ^ ¥'(") = C2 (49) 

Since the set of differential equations for the geodesic line is invariant with respect to the translation 
for variables t, z and (f, we can assume that t{Q) = 0, (p{0) = and z{0) = 0. We choose the 
other initial values for the geodesic lines randomly. For example we can choose the initial values 
r(0) — 1.33, r(0) = 0.78, t(0) = 0.81 and (j){0) — 0.56. From these values we have the exact numbers 
/(0, 1), / (l, —j), Cl and C2 from Equations (|46ll49p . We use the MAPLE software, dsolve package, 
dverklS method to solve numerically the geodesic equations with these initial values on the interval 
[0,2]. We set the absolute and relative error for the solution to be 10^*. We compute the quantities 
of the left hand side in the Equations ()46M9p at the points u — k ■ 0.005 (1 < fc < 400), and we 
compare the computed values to the exact ones. In this case the biggest error for Equation is 
El = 0.5 X 10"^ for Equation g?]) is £^2 = 0.5 x 10"^ for Equation gH]) is £^3 = 0.4 x 10""^ and for 
Equation ([i^ is £4 = 0.54 x 10~^. The results of 20 simulations are shown in Table [TJ 

It is clear from the table, that the errors are acceptable with respect to the computational precision. 
This gives a numerical verification of the Theorem. 

2.4 Some simple consequences 

If we choose the parameters in the matrix M to be (pi, P2) = (2, 0) and (pi, P2) = (0, 1) then the first 
components of the Equation are 



sinhV) . ^ /(2,0) 
V2 cosh^(r) 



(50) 



t + V2sinh2(r)(^ = /(0, 1). (51) 
Solving these Equations for t and (p we have 

/(0,1) (52) 



2/(2,0) 



cosh^ (r) 



sinh^(r) V cosh^(r) 

So the functions t and ip depend just on the initial conditions and on r. Moreover we can define a 
critical radius Rt — archy^^^^ (if exists) from the initial conditions, where i — Q. \i r > Rt then 
the particle moving backward in time and if r < Rt then moving forward. 
This leads us to the following Lemma. 
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Table 1: 



r[{)) 


r({)) 


t[U) 


0(Oj 


-C/i X lU 


11/2 X iU 


TP \/ 1 o8 


TP \/ 1 n8 

i!/4 X iU 


L.OO 


U. / o 


U.oi 


U.OO 


U.O 


U.O 


A 
4 


O 

U.04 


11/1 


n 01 


U./D 


1 QV 


7 

u. * 


U.D 


Z.D 


O /lO 
U.4Z 


n K 
U.O 


— u.y ( 


U.D < 


n 1 /I 
U.14 


U.zo 


U.io 


O 7Q 


q 
o 


U.Do 


1 


n 07 
U.y < 


— U.o4 


U.z 


U.io 


A /I 
0.4 


4 


n 1 

U.i 


1 /I Q 


1 00 


— U.O / 


n A 
U.4 


n 1 o 
U. iy 


O Qfi 
U.OO 


Q 
U.O 


U.o4 


— U.OD 


l.oD 


111 
i.li 


U.z 


n o 
U.z 


o c 
U.O 


n oo 
U.zo 


U.Oz 


u.yo 


U.O 


— u.oy 


O 07 

U.U / 


U.U4 


OA 1 Q 
ZD. iO 


A 





n 1 Q 
U.io 


U. / 


U.O / 


U.O 


O Q 

U.O 


z 


o o 
u.z 




U.oz 


1 Ofi 
l.ZD 


n A7 
U.D / 


n 1 
U.i 


u.i 


O 1 Pi 
U.io 


O OA 
U.UO 


nil 

U.ii 


— U.oz 


1 


n OQ 
— U.Uo 


1 
U.i 


1 

u.i 


O 1 
U. i 


Oil 
u.ii 


1 OQ 


n Qi 


U.D ( 


— U.Oo 


n /I 
u.4 


U.ZD 


1 7 
i. / 


O A 
U.D 


n 1 /I 


— U.oZ 


1 A 


n f^i 
—u.oi 


n 1 

U. i 


n o 
u./ 


O 1 7 
U.i / 


O 1 O 

u. iy 


U.OO 


u.yz 


u.uz 


U A*? 

U.'iO 


U.U4 


n no 


o.oz 


1 

1 


1.24 


0.79 


0.69 


1.25 


0.8 


0.8 


5 


0.8 


0.15 


-0.17 


1.49 


0.14 


0.1 


0.1 


0.15 


0.3 


0.3 


0.45 


0.35 


0.97 


0.09 


0.06 


0.11 


0.28 


0.41 


0.69 


1.49 


-0.08 


0.1 


0.1 


0.4 


0.19 


0.7 


-0.2 


1.41 


1.48 


0.1 


0.2 


0.06 


0.04 


0.2 


0.17 


1.14 


-1.00 


0.1 


0.1 


0.12 


0.04 


0.36 


-0.84 


1.31 


-1.08 


0.2 


0.1 


0.3 


0.18 



Lemma 2.4. ///(2,0) > and /(0, 1) < then i{u) > for every u e Dom7. If f {2,0) < 
/(0, 1) > then t{u) < for every u G Dom7. 

Let us compute the inverse of the matrix M 

(o A ui'o \ crcos((p) tanh(r) a sin((^) tanh(r) \ 



1 



8-4cosh(2r) 


8\/2 




V 



V2 
K sin((^) 

2cr cos((/3) 
sinh(2r) 




V2 

K C0s((/7) 

— 2crsin((^) 
sinh(2r) 











where 



K = pi + (pi + 4^2) cosh(2r) and a = pi + 8p2 + Pi cosh(2r). 
From these Equations proves the following Lemma. 
Lemma 2.5. We have correspondences betiueen the functions r and (f 
r = ci sin((^) + C2 cos((p) 



sinh(2r) 

Which means that the function f is bounded 



(ci cos(<p) - C2 sin((p)) + /(2, 0) 



2\/2 
cosh^ (r) 



r < 



(54) 



(55) 

(56) 
(57) 

(58) 
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and we have a bound for tp 



< 



sinh(2r) 



+ /(2,0) 



2^/2 



(59) 




From Equations ()53|56p we get 



dr 



ci sin(i^) + C2 cos{ip) 



(60) 



d if 




) 



After integration we have the foUowing Lemma. 



Lemma 2.6. There is a real parameter m such that 



C2 sin((p) — ci cos((^) 



2V2 



(/(2, 0) cosh(2r) - /(0, 1) cosh^{r)) + m 



(61) 



sinh(2r) 



holds. 
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